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Abstract. We analyze the quantum dynamics of the time-dependent elliptical 
billiard using the example of a certain breathing mode. A numerical method for 
the time-propagation of an arbitrary initial state is developed, based on a series of 
transformations thereby removing the time-dependence of the boundary conditions. 
The time-evolution of the energies of different initial states is studied. The maximal and 
minimal energy that is reached during the time-evolution shows a series of resonances 
as a function of the applied driving frequency. At these resonances, higher (or lower) 
lying states arc periodically populated, leading to the observed change in energy. The 
resonances occur when the driving frequency or a multiple of it matches exactly the 
mean energetic difference between the two involved states. This picture is confirmed 
by a few-level Rabi-likc model with periodic couplings, reproducing the key results of 
our numerical study. 
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1. Introduction 

Recently, the dynamical properties of classical time-dependent two-dimensional billiards 
have attracted major attention, especially in the context of Fermi acceleration 
[HEllHlllllHlinilTllH], which is defined as the unbounded energy gain of an 
ensemble of particles exposed to some external driving force [9]. Concerning the 
quantum dynamics of time-dependent billiards, there exist several studies investigating 
the quantum version of the one-dimensional Fermi-Ulam model (or variants of it) 
[IOl[IIl[I3|13l[Tll[T5l[T6l[I71[T8l[l9l|2^ but, to our knowledge, only two studies 
[221 123] investigate time-dependent billiards in higher dimensions. 

Since already the one-dimensional Fermi-Ulam model involves time-dependent 
Dirichlet boundary conditions which are difficult to treat, both analytically and 
numerically, most of the works analyze which (non-periodic) movements of the wall 
allow for exact solutions: Linearly expanding or contracting wall motion is considered 
in this context in Refs. and [11]. The authors of Refs. [121 US] find that exact 
solutions are not only possible for a linear wall motion, but also for a time-law of 
the form l{t) = V af^ + 2bt + c (with appropriate real constants a,b and c). The same 
time-law, that allows for analytic solutions, is found from a different perspective in Refs. 
[T5|[T8] and by means of a supersymmetry formalism in Ref. [20]. If additionally specific 
electro-magnetic fields are superimposed, exact solutions can be obtained for arbitrary 
time-dependencies of the wall motion [19] . Periodic driving laws are considered in Refs. 
[TT] and [21]. The Floquet or quasienergy spectrum is found to be pure point like for most 
laws of the wall oscillation, resulting in a recurrent dynamics and thus to bounded energy 
growth [11]. Only certain (non-smooth) driving laws yield a continuous quasienergy 
spectrum and thus allow for an unbounded acceleration, for example via resonance 
excitations, similarly to the quantum kicked rotator [23] . Recently, the authors of Ref. 
[2T] studied the quantum FUM numerically by expanding the wave function in terms 
of the instantaneous eigenstates of the corresponding static system. They find that the 
dynamics of the time-dependent expansion coefficients is irregular in an intermediate 
frequency regime, whereas for very low and very high frequencies a periodic behavior of 
the expansion coefficients is obtained. 

Concerning higher dimensional billiards, the one-pulse response of a two- 
dimensional stadium billiard to a deformation of the boundary has been studied [23] 
by analyzing the evolving energy distribution. For small wall velocities, the spreading 
mechanism of this distribution is dominated by transitions between neighboring levels, 
while this is not the case for non-adiabatic wall velocities. In Ref. [22], the radially 
vibrating three-dimensional spherical billiard is investigated. The authors claim [|| that 

I Among others, their arguments are based on the correspondence between the classical and quantum 
3D driven spherical billiard. By doing so, they state that in the classical system, the angular momentum, 
which is a constant of the motion in the static case, gets destroyed by the driving. This is not correct, 
since it is known, see for example Ref. [2], that in the radially oscillating two-dimensional circular 
billiard the angular momentum is still conserved and the arguments can be generalized straightforwardly 
to the 3D spherical billiard. 
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only superpositions of two or more states that share the same common rotational 
symmetry yield quantum chaos, since the orthogonality relations of the instantaneous 
eigenstates allow in any other case a reduction to a time-independent one degree of 
freedom Hamiltonian which cannot be chaotic. 

In the present work, we investigate the quantum version of one of the most 
studied classical two-dimensional time-dependent billiards, the elliptical billiard with 
harmonically oscillating boundaries. Already the classical system shows several 
astonishing results: the driven elliptical billiard exhibits (tunable) Fermi acceleration [6] 
and the corresponding diffusion in momentum space experiences a dynamical crossover 
from sub- to normal diffusion [251 ESj. A more technical appealing aspect of using the 
elliptical billiard is that the quantum version of the static billiard can be solved exactly 
in terms of the Mathieu functions, which allows an intuitive analysis in terms of the 
instantaneous eigenstates when driving the system. In particular, we will study the 
time-evolution of the energy in case of an initial eigenstate of the corresponding static 
billiard. The maximal and minimal energy that is reached as a function of the driving 
frequency show a series of resonances which we explain in detail by means of a detailed 
population analysis. We develop a Rabi-like model with periodic couplings which nicely 
reproduces the results of our full numerical simulations, confirming the intuitive picture 
of the level dynamics of the driven elliptical billiard. 

The work is structured as follows. In section |2] we summarize how the eigenstates 
of the static elliptical billiard can be obtained by using the Mathieu functions and 
discuss the symmetry properties of the eigenstates. Subsequently, we tackle the problem 
of how to solve the driven elliptical billiard numerically in section [31 The results of 
the simulations are presented in section H] and interpreted, explained and modeled in 
section |5l Finally, a short summary and outlook is provided in section |6l 

2. Static elliptical billiard 

As a precursor for the analysis of the time-dependent elliptical billiard, we briefly study 
the static elliptical billiard. In particular, we will use in the analysis of the time- 
dependent system projections on eigenstates of the static ellipse and the symmetry 
properties of these eigenstates. To find the quantum mechanical eigenstates \E'„(a;, 
of a static elliptical bilhard, we have to solve the corresponding stationary Schrodinger 
equation 



where E is the eigenenergy belonging to \&(a;, y) and the Hamiltonian H{x,y) is given 




(1) 



by 




(2) 
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where n denotes the mass. The potential V{x,y) is zero inside and infinite outside the 
static eUiptical bilhard 

'<i-$fi (3) 

[ OO lf^ + ^ > 1. 

This means the eigenvalue equation ([1]) has to be solved under Dirichlet boundary- 
conditions 

f 2 2 ^ 

^{x,y) = for l{x,yy eM'\^ + ^ = l\, (4) 

where a and b are the semi-major and semi-minor axis of the ellipse, respectively. 

Due to the boundary condition of Eq. @, solving the eigenvalue equation ([1]) is 
actually not trivial and involves some subtleties [271 1211 [291 EOl EH |32l [33]. Here, we 
follow the discussion provided in Ref. [32] and additionally use results presented in Refs. 
[32] and [31]. We start by transforming from Cartesian to elliptic coordinates: 



X = / cosh ^ cos ?7, < ^ < ^0 (5 a) 

y = /sinh,^sin?7, < rj < 2tt, (56) 

where / = \/a^ — 6^ is the semi-focal distance of the ellipse and = arctanh(6/a). The 
stationary Schrodinger equation reads 

H{^,r^M^,r^) = Ecpi^,v), (6) 

where (p{C,,'ri) is an eigenfunction with energy E and H{^,r]) is the Hamiltonian: 

H{^,v) = -^^i^,v) + V{^,v) (7) 
The potential energy V{C,,ri) in elliptic coordinates is now 



if^ < ^0 
OO if^ > ^0, 



As a consequence, the eigenfunctions (p{C,,f]) have to satisfy the Dirichlet boundary 
condition (p{^o,r]) = and also the periodicity condition ip{^,ri) = (p{^,r] + 27i). The 
Laplacian operator in elliptic coordinates is given by 

2 / 92 ^2 \ 

^ (^' ^) = j2 (cosh 2^ - cos2ri) [oe ^ W ) ' 
The two-dimensional (2D) Schrodinger equation inside the elliptical billiard (i.e. ^ < ^0) 
then simplifies to 



g2 Q2 f^2 p 

de^W^ ^(cosh2e - cos2r/) 



¥>{C,V) = 0, (9) 
with k"^ = 2fiE/h?. Obviously, ip{C,,ri) can be separated with the standard ansatz 

^i^,v) = Rmiv) (10) 
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into a radial (-R(O) ^ind an angular (6(77)) part. Plugging this separation ansatz in the 
Schrodinger equation ([9]), we obtain two ordinary differential equations: 

R"{0 - (« - 2gcosh20/2(0 = 0, (11a) 
e"{ri) + (a- 2q cos 2ri)e{ri) =0, (115) 

where a is the separation constant and q is the dimensionless, rescaled energy 

Equations (lllap . flll6p are the standard form of the Mathieu equations, Eq. fillip is 
the ordinary Mathieu equation (OME) and Eq. f ilial) is the modified Mathieu equation 
(MME). The corresponding solutions are the ordinary and modified Mathieu functions 
(OMF, MMF), respectively. Note that the change of variables rj = iC, transforms the 
OME into the MME. Even though the ansatz ffTOj) decouples the Schrodinger equation 
([9]), yielding two ordinary differential equations, the separation constants a and q do 
not decouple [3l], as we will see in the following. 

The physical solutions of the OME have to be periodic, i.e. they can be expanded in 
a Fourier series (actually, Floquet theory guarantees the existence of periodic solutions 

m) 

edr]) = I " ^^=0 ^°^(^^) (13) 

' \ sez(r/, q) = Bk{l, q) sm{krj) 

For a fixed q, it is known that only certain values of a allow for periodic solutions. 
These eigenvalues are called the characteristic values ai and f3i, where / is the order 
(/ > for cei{ri) and / > 1 for sei{r])). Since the Mathieu equation is of Sturm-Liouville 
type, all eigenvalues are real, positive and can be ordered ao < /3i < cti < /32 • • •• In 
the q X a plane, the functions ai{q) and Pi{q) are curves that do not intersect. The 
cei(?7, q) and se;(?7, q) with even order have a period of vr, the ones with odd order have 
a period of 27r. The order / specifies the number of zeros of ce;(?7, q) and se;(?7, q) in the 
interval 77 G [0,7r]. The expansion coefficients Ak{l,q) and Bk{l,q) are determined by 
the recurrence relations for the Mathieu equations |35j . 

The solutions of the MME Eq. ( lllap . i.e. the radial part, can be obtained from the 
solutions of the OME by the mentioned change of variable rj = iC,, yielding 

[ Sei{^,q) = ZT=iBk{l,q)smh{kO. ^ ' 

Note that the expansion coefficients Ak{l,q) and Bk{l,q) are the same as the ones in 
the solutions (fT3l) of the OME. The solutions of the MME can also be viewed as an 
independent eigenvalue problem for the g's (with fixed a and /3), now Ri{^) has to the 
satisfy the Dirichlet boundary condition Ri{^o) = 0. 

The eigenvalues associated to a and (3 are denoted by g^(a) and q°.{P)^ where the 
index o, e refers to the even and odd MMF and r > 1 is the rth zero of the MMF, not 
counting a possible root at = 0. In the g x a plane, the functions q^{a) and q°.{P) are 
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curves that again do not intersect. The solution f{C,,f]) of the Schrodinger equation 
is the product ffTOl) of the OMF and MMF. Both eigenvalue problems have to be solved 
simultaneously, since a (or /3) and q'^ (or q°) in the OMF and MMF of the solution 
^{i,rj) cannot be chosen independently. The pairs {a,q) (or (/3,g)) that satisfy both 
eigenvalue problems simultaneously are of course the crossing points of the family of 
curves ai{q) and (3i{q) with the q^ia) and q°{/3). Thus, there are two types of crossing 
points that correspond to solutions of the simultaneous eigenvalues problem: 

(i) ar{q) with q^i^:) (even modes) 

(ii) f3r{q) with q^if^) (odd modes) 

These two types of crossing points are related to the symmetries of the solution ip{C,, rj). 
This yields the eigenf unctions ipi^riC^v) 

^lA^, V) = Cei{^, qlr)cei{7], gf j (15a) 

where qf^ = qr{c(i) and correspondingly = q°{(3i). The indices / and r are the angular 
and radial quantum numbers of the even and odd modes, respectively. / is the order of 
the OMF and r specifies the number of zeros of the MMF in the interval ^ G [0, ^o], not 
counting the zeros at the origin i.e. gf^° is the r-th zero of the MMF of order /. Since 
the MMF has to be zero at it follows that r is greater or equal to one. According 
to Eq. f|T2|) . the separation constant q is directly proportional to the energy and the 
energy eigenvalues El'° corresponding to the eigenfunctions ipl'°{^i v) can be written as 

7-ie,o 2/l gQ 2h gQ 

where e is the eccentricity of the ellipse. 



Symmetries 

The symmetries of the eigenstates '^1'°{x,y) of the elliptical billiard are determined by 
the OMF Qiijf) (IT^ [21]. For an even OMF, i.e. for ceiijf), the wavefunction \Ef^j,(x, y) 
is symmetric with respect to reflections at the x-axis (\E'^,,(x, y) = ^^^,,(0;, — y)), this 
symmetry is denoted by tTj^ = +1. Odd OMF sez+i(?7) imply an antisymmetric eigenstate 
^°r(^'2/) with respect to the x-axis {^^^{x,y) = — \E'^^(a;, — y)), denoted by iTy = — 1. 
The symmetry of ^^'°(x, y) with respect to the y-axis is fixed by the angular quantum 
number /, which is the order of the OMF. For even Mathieu functions, \E'fj,(x,?/) is 
symmetric with respect to reflections at the y-axis {'^°^{x,y) = '^°^{—x,y)) if / is even 
and antisymmetric (^'°^(a;, |/) = — \l/°^(— x, y)) if / is odd. Naturally, this symmetry is 
denoted by tTx- For odd Mathieu functions, \E'°,,(a;, |/) is symmetric to the y-axis if / 
is odd and antisymmetric if / is even. Overall, this yields four possible combinations 
concerning the parity of the eigenstates \E';^'°(a;, y), which are summarized in the following 
hst: 

(i) Tiy = +1, Tlx = +1 ^ r(^' y) — r(~^' Gveu solutious with even / 
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Table 1. All eigenstates 'i>1'°{x,y) of the elliptical billiard {a ^ I, b — \/0.51 => 
eccentricity e = 0.7) with E < 50 (in units of /i-l). The eigenstates are ordered 
according to the energy. TTy (tt^^) denotes the symmetry of ^^'°(a;, y) with respect to 
reflections at the x-axis (y-axis). 



Label 




/ 


r 


even/odd 


■Ky 




1 


4,267 





1 


e 




+ 


2 


9,058 


1 


1 


e 




— 


3 


12,577 


1 


1 


o 


— 


+ 


4 


15,993 


2 


1 


e 




+ 


5 


19,358 


2 


1 


o 






6 


25,061 


3 


1 


e 






7 


25,895 





2 


e 




+ 


8 


27,998 


3 


1 


o 




+ 


9 


35,156 


1 


2 


e 






10 


36,178 


4 


1 


e 




+ 


11 


38,511 


4 


1 


o 






12 


44,040 


1 


2 


o 




+ 


13 


46,406 


2 


2 


e 




+ 


14 


49,199 


5 


1 


e 







(ii) Tiy = +1, Tlx = —1 ^ ^Fr(^' y) — ~^Fr(~^' Gven solutions with odd / 

(iii) TTy = —1, TT^. = +1 ^ \E'°^(x, y) = — ^°^(— x, —y) odd solutions with odd / 

(iv) Tiy = —1, Tlx = — 1 ^ = X, —y) odd solutions with even / 

The four parity combinations correspond to the characterization of the symmetry 
reduced quarter elliptic billiard. However, the boundary conditions along the coordinate 
axes of the quarter billiard have to be adjusted according to the different parities: 
Dirichlet boundary conditions are required along the x-axis (y-axis) for Tiy = —1 
[tTx = —1) and Neumann boundary conditions along the x-axis (y-axis) for Tiy = +1 
(vr, = +1). 

3. Driven elliptical billiard 

To study the time evolution of an initial state '^^{x^y) := '^{x,y,t = 0) in 
the harmonically oscillating elliptical billiard, we have to solve the time-dependent 
Schrodinger equation 

^^ d^{x^y,t) ^ ^^^^ ^^^^^^ ^^^^ 

where the Hamiltonian H{x,y,t) is given by 
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and the potential V{x, y, t) is zero inside and infinity outside the time-dependent 
elhptical bilhard: 



„2 



if ^ + t^It-y < 1 

I o« + b%) > 1- 

This means the Schrodinger equation (|T7I) has to be solved under Dirichlet boundary 
conditions 

vl/(x,2/,t) = for |(x,2/)^gM2|^ + ^ = 1, Vtj, (20) 
where a{t) and 6(t) are determined by the driving law of the elliptical billiard. 

Transformations 

The explicit time- dependence of the boundary conditions is very difficult to be treated by 
standard numerical techniques. Therefore, to remove the explicit time-dependence of the 
boundary conditions, we apply the following time-dependent coordinate transformation, 
whose ID variant has been successfully used to remove time-dependent boundary 
conditions [ini[12l[l3l[Tll[l8l[l9l[20l|2l]inlD systems 

yielding the time-dependent Schrodinger equation 

^^ d^mm.t) ^ Hm,at),tMvit),m,t). (22) 

The advantage of applying this time-dependent coordinate transformation is that the 
boundary condition in the new coordinates 77, ^ is extremely simple: the wave function 
\E'(?7,,^,t) has to vanish on the circle given by 77^ + .^^ = 1. The prize we have to pay 
is that now the coordinates themselves are explicitly time- dependent, rj = rjit) and 
^ = ^(t). This has to be taken into account when applying the differential operator 
d/dt on '${T]{t),^{t),t) in Eq. (122|) . The additional terms resulting from the left hand 
side of Eq. f l22p can be put into an effective Hamiltonian, yielding 

^' ^' ~ 2/i \a^{t) dr]^ ^ b^{t) de 
and the corresponding time- dependent Schrodinger equation 



ih- 



dt 



= H%7^,^,t)^{r],^,t). (24) 

ri,(^=const. 



To remove the terms proportional to ihri-^ and ih^^, we apply the following time- 
dependent unitary transformation, yielding the new wave function A(?7, ^, t) 

^{r],U) = ^{v,t)-r{U)-Hv,U) (25) 

2 ( / • 2 I i,i,f-2^ 



vl>(r/,e,t) = -^exp( ^(aar/^ + &&e')) ■A(r/,e,t). (26) 
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With this transformation we obtain the Schrodinger equation 

^^9A(|e^ = w{v, e, t)A{v, e, t), (27) 

with the effective Hamiltonian 

•H^(r/,e,t) = -— ' 



2/i \a^{t) 9r/2 b^{t) de 
+^/ia(t)a(t)r/2 + ]^^^h{t)h{t)e (28) 

and the time-independent boundary condition A{r],C,,t) = for 77^ + = 1. This 
effective Hamihonian can be interpreted as a two-dimensional anisotropic harmonic 
oscillator with time-dependent frequencies and time-dependent masses and the above 
boundary condition. 

We expand A(?7,^,t) in terms of the eigenfunctions $n,m. of the (unit) static 
circular billiard, since this ansatz automatically fulfills the boundary condition. The 
eigenfunctions $n,m of the static circular billiard are best described in polar coordinates, 
so we transform rj and ^ to polar coordinates: 

7^ = rcos0 ^ = rsin0 ^ A(r = 1, 0, t) = 0, V0, t (29) 

The Hamiltonian (l28l) in polar coordinates reads 

/ . 2 i 2 sin cos d 



2na^{t) \ dr"^ dcj) 

2 sin cos (9^ cos^ d cos^ d"^ 
+ -^r^ + 



2 /)^2 



r 9r90 r dr r 

/ 2 _|_ 2 sin cos (j) d 2 sin cos d"^ 
2/x&2(t) V Sr^ d(j) r drdcj) 

sin^ 4) d sin2 9^ \ 1 2 2 



-fia{t)d{t)r cos 



r 9r 90^/ 2 

+ ^/i6(t)6(t)r2sin2 (30) 

and the corresponding Schrodinger equation is 

^^ gA(^0,t) ^ ^^^^^^ _ ^g^^ 

Ansatz 

As already written, we expand the wave function in terms of the eigenfunctions 
^n,m{'^, 0) of the static circular billiard. To carry out numerical simulations, this 
expansion has to be truncated, yielding 

N M 

A(r,0,t) = ^ J2 c„,„(t)$„,„,(r,0), (32) 
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with time-dependent expansion coefficients Cn^mit), where n is the radial and m the 
azimuthal quantum number. The time- independent eigenf unctions of the 

static circular billiard factorize 

<^n,m{rA) = Rn,m{r)-QM)- (33) 

The normalized radial and azimuthal (angular) functions are given by [36l 



V2 

Rn,m{j') = 'J 77 ^Jm{km,n'^) (34) 
'J m+l\'^n,m ) 

QM) =^e'""^ (35) 
where Jm is the Bessel function of the ffist kind of order m and km,n is the nth zero of 

To determine the time-dependent expansion coefficients Cn,m{t), we insert the ansatz 
( l32l) into the time- dependent Schrodinger equation ( |3T]) . project with the bra (A| 

{A\ih^\A) = {A\H^{t)\A) (36) 

and thus obtain a set of coupled ordinary differential equations (ODE) for the Cn^m{t)- 
This coupled ODE system can be written in the following form, for details of the 
derivation see Ref. [38] 

iflCn,m{t) = 5]]{c„',„(t)4mn'(^) 

n' 

+C„',^-2(t)41n'(^) + Cn',™+2(t)41„'(^)}- (37) 

The d^^^^,{t) can be written in the following form (exemplarily shown for d^^^^,{t)) 

where the g^^\t) are simple time-dependent functions of a{t) and 6(t), whereas the 

fnmn' time-indepcndent. For a definition of the (i^*^„/(t), the f^mn' ^^^1 the g'^^^t), 

CO 



see [Appendix A[ More specific, the are composed of integrals over products of 

Bessel functions, which in general cannot be evaluated analytically. However, since 
the corresponding integrals are time-independent, they have to be evaluated only once, 
which can be done numerically with high precision. Thus, by introducing a linear index 
(n, m) -v^ i the coupled ODE system ( !37l) is in the canonical form and can be solved 
numerically by standard techniques. 



Observables 

The time-evolution of the expectation value of the energy of a certain initial state 
propagating in the driven elliptical is given by 

E{t) = {m\H{t)\m)- (38) 

To calculate E{t) at a given time to of the original laboratory system, we cannot directly 
use the wavefunction A(r, 0, tg) fl32|) together with the Hamiltonian 'H'^(r, 0, to) fl30|) . 
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Instead we have to apply the following steps: First, the inverse of the transformations 
fl2T]) . fl25|) and (129|) operate on A(r, to) and we obtain 



xij(^x,y,to) = ^— e^('^l'=*o'''/''^*°)+^l'='o^'/*^*°^) 
yob 



J2cn,m{to)^n,m ( + ^rctan H 

n,m \ / 



(39) 



Secondl, we apply the coordinate transformation as given in Eq. (12T|) . but now as a 
time- independent coordinate transformation, which just depends parametrically on the 
time to'. 

a{to) o{to) 
The wavefunction and the Hamiltonian are then given by 

'^Cn,^ito)^n,m ( ^/ + , aictau (-) ) , (41) 

n,m ^ \^ / / 

Finally, changing again to polar coordinates yields 

vi/(r,0,to) = e^^'('^'^+(^^-'^'^)^^'^'<^) 5^c„,^(to)$n,m(r,0), (43) 

n,m 

and the corresponding Hamiltonian 'H(r, t) is given by Eq. (I30l) if we set a(t) = 6(t) = 
0, i.e. 

H{r,ct^,t)=n%r,c^,t)\.,^;^,. (44) 

Another quantity of interest is the spectral composition of the energy in terms of the 
population of the instantaneous eigenstates of the elliptical billiard. By 'instantaneous 
eigenstates' we mean the following: Consider the boundary of the elliptical billiard for 
a certain time instant. If we now take this particular billiard configuration, treat it as a 
static ellipse and calculate the corresponding eigenstates assuming Dirichlet boundary 
conditions, we obtain the instantaneous eigenstates. If we denote the eigenstates by 
(\&j(t) I {i labels the degree of energetic excitation), the population of an instantaneous 
eigenstate is then of course simply the projection of the wave function '^{x,y,t) onto 
this state, i.e. pi{t) = \ (^'i(t)|^'(t)) |. 

It is illustrative to analyze how the energy of the instantaneous eigenstates changes 
during one oscillation of the elliptical billiard. From Eq. (fT6l) it can be seen that the 
energy of the eigenstates of the elliptical billiard depends (for a given semi-major axis a) 
on the eccentricity e, so Ei = Ei{e). Since the periodic (with period T = 2'k/u) driving 
law a(t), h{t) changes the eccentricity, we may write for the energy of the instantaneous 
eigenstates Ei{(), where ( is the phase of the oscillation, measured in terms of T, i.e. 
( e [0, 1]. Ei{() for energetically lowest 14 states is shown in figure [H There are several 
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Figure 1. (Color online) Energy -E'i(C) (in units of fi^/ji, the units will be omitted 
in the following for simplicity) of the instantaneous eigenstates as a function of the 
phase C of the elliptical billiard. The solid lines represent the energetically lowest five 
eigenstates with symmetries tt^ ^ TTy = +1. The inset shows the crossings of the states 
six and seven. 




10 



20 30 
time t/T 



1: l=0,r=1,e 
-4: l=2,r=1,e 



40 



50 



Figure 2. (Color online) Evolution of the energy E{t) for a driving frequency of a; = 5, 
when starting in the instantaneous ground state (first state) and in the fourth state. 



crossings, note that these are real crossings and not avoided crossings, see the inset. 
Such crossings are possible since the static elliptical billiard is an integrable system, 
thus the corresponding eigenstates posses 'good' quantum numbers, in agreement with 
the theorem of Wigner and von Neumann [39] . 
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Figure 3. (Color online) Time-evolution of the projections Pi{t) onto the 
instantaneous eigenstates when starting in the fourth state, i.e. P4(t) = 1 with a driving 
frequency of w = 5, of. figure [2j Besides p4, now gets considerably populated, note 
the large beating period T;, > T of both coefficients. 



4. Results 

In the following, we will take as an initial state \E^o = ^(^ = 0) an eigenstate of the static 
elliptical billiard, i.e. \l/o = ^^"(a^,?/), and let it evolve in the driven ellipse. We choose 
for the driving law of the elliptical billiard a{t) = ao + csin{ujt) and b{t) = 6o + csin(ajt). 
For the equilibrium positions of the semi-major and semi-minor axes we use ao = 1, 
Bq = a/0.51, resulting in an eccentricity of eo = 0.7. The driving amplitude c is set to 
0.1, whereas the frequency u will be varied. 

In figure [21 the evolution of the energy E(t) is shown for a driving frequency of 
u = 5 and two different initial states, namely the ground state (first state) and the 
third excited state (fourth state), i.e. \E'o = "^oii^^u) = ^21(^52/)' respectively. 

These are the first two states that have the same symmetry, namely vr^ = tt^ = +1, 
see table [H When starting in the instantaneous ground state (at t = 0), the wave 
function \Ef(t) follows the motion of the elliptic boundary adiabatically. In other words, 
the projection pi{t) of \E'(t) on the instantaneous ground state is approximately one, 
Pi(^) ~ I (^0 I ~ 1- The periodic oscillations of the energy E{t) are then 

simply present due to the fact that the energy of the instantaneous ground state changes 
periodically (approximately sinusoidal), cf. figure [1] 

The situation is different when starting in the fourth state (\l/o = ^'2 li^yV))- Even 
though E{t) performs sinusoidal- like oscillations corresponding to the energy oscillations 
of the fourth instantaneous eigenstate, there is a small modulation (beating) on top, 
with a period Tb > T = 2it/u) much larger then the period of the applied driving law. 
This large-periodic modulation is better visible in the time-evolution of the population 
coefficients Pi(t), which are shown in figure [3l Since \E'o = ^21(^5^)5 P4(0) = 1 and 
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Figure 4. (Color online) Evolution E{t) for a driving frequency uj ~ 10, when starting 
in the instantaneous ground state (first state) and in the fourth state. Whereas the 
fourth state shows a pronounced beating phenomenon, this is not the case for in the 
ground state. 
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Figure 5. (Color online) Evolution E{t) for a driving frequency of w = 12, when 
starting in the instantaneous ground state (first state) and in the fourth state. Both 
curves show a pronounced beating phenomenon. 

all the other pi{0) are zero. The fourth instantaneous eigenstate stays populated 
(p4(t) > 0.98, Vt) throughout the whole time-evolution, however - besides small, fast 
oscillations with the period T of the applied driving law - there is now a beating with a 
much larger period Tf,. Correspondingly, other states get noticeably populated, primarily 
Pjit) and Pi{t) are now greater than zero, with pj having the same large beating period 
Tb. Note that the instantaneous eigenstates corresponding to pi, p4 and pi share the 
same symmetry properties tTj. = tt^, = +1, cf. Tabled! 

When choosing other low energy eigenfunctions as an initial state, e.g. \E'o = 
'^l i(x,y), 'i/° i{x,y), \E'2,i(^5^) = ^3,1(2^)^)! the result is very much like the one 

when starting with \E'o = ^0 li^^ u)-, in the sense that the energy E{t) performs periodic 
oscillations matching the sinusoidal-like oscillations of the corresponding instantaneous 
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eigenstate. This behavior is also observed for smaller driving frequencies < 4), now 
also when starting in the fourth state, i.e. \E'o = ^2 v)- 

For larger driving frequencies, the beating phenomenon becomes even more 
pronounced. In figure HI a driving frequency of w = 10 is applied to \i/o = ^0,1(2^) J/) 
and \l/o = ^2i(^5l/)- For \E'o = "^oii^yV) resulting E{t) shows besides the large 
amplitude beating some irregular fluctuations, however E{t) stays tightly within a well- 
defined band. Now for \&o = ^21 (^'2/) beating is much more pronounced. Again, 
the period Tb of the beating is much larger than the driving period T, note however 
that it is not the same beating period we obtained in the case of a; = 5, cf. figure [31 
Like in the case u = 5, the time-evolution of the population coefficients, see figure El 
reflects the beating phenomenon. Again (for \&o = ^2 1(^5 2/))' ^^e coefficient is the 
one right behind p^ that gets populated the most, note that at some times even pi > p^. 
Other states with noticeable population are pio, pi and pi3, again these are all states 
with symmetries vTa. = vr^ = +1, cf. Table [H 

By increasing the driving frequency further to a; = 12, the beating phenomenon is 
now observed equally when starting in the instantaneous ground state (^'o = ^0 u)) 
and when starting in the fourth state (^'o = ^2 1(^5 1/))' see figure [5l Unlike in the 
case u = 10, the large-periodic modulation (in the case \E'o = "^21(^^11)) considerably 
lowers rather than rises the energy E(t). Consequently, the state besides p4 that gets 
substantially populated is now pi and not the higher energetic p^, see figure [71 Even 
more, pi reaches almost one, i.e. almost the whole population is periodically transferred 
between p4 to pi. 

The beating phenomenon gets destroyed for sufficiently high driving frequencies, 
see figure[8l where E{t) for u = 100 with \E'o = "^o^iix, y) (pi(0) = 1) and \E'o = ^l^i{x, y) 
(p4(0) = 1) is shown. Now a considerable number of states with the same symmetry 
as the initial state get noticeably populated. Consequently, the energy does not stay 
in a narrow interval as in the case of small driving frequencies, but rather fiuctuates 
irregularly over a wide range. More specific, E{t) increases rapidly within the first 
few oscillations of the elliptical billiard, then slowly grows further until t ^ 60T and 
saturates for t > 60 T. This means that the energy stays bounded, unlike in the 
corresponding classical system [25l [26], there is no Fermi acceleration in the quantum 
version. This is due to dynamical localization: the wave function gets exponentially 
localized in momentum space, thereby stopping the energy diffusion. 

In general, when starting at t = with an instantaneous eigenstate with energy 
Ei, the maximal -Emax = max(£'(t)) and minimal energy -Emin = min(£'(t)) of \E'(t) is a 
good indicator whether the time-evolution of "^{t) is adiabatic or not. If the evolution 
of \i/(t) is adiabatic, -Emax.min are simply given by the corresponding eigenvalues of the 
instantaneous eigenstate, more specific, in this case {u = 1) -Emax = Ei{( = 3/4) and 
-Emin = Ei{( = 1/4), cf. figure [H As soon as the time-evolution is not purely adiabatic, 
the beating phenomenon sets in and either -Emax is raised, or -Emm is lowered, or both. 
This can be seen in figure [HI where -Emax and -Emin are shown as a function of the 
driving frequency u when starting with the fourth state (\I'o = ^01(^5 1/))- -^o^ small 
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time t/T 

Figure 6. (Color online) Time-evolution of the coefficients Pi{t) corresponding to the 
upper curve (state 4) of figure |4] (cj = 10). Again p4 and pr dominate, however, now 
further states {pio,Pi,Pi3) contribute. 




time t/T 

Figure 7. (Color online) Time-evolution of the coefficients pi{t) corresponding to 
figure[5](u; = 12), when starting in the fourth state. Unlike in the case uj = 10, besides 
Pi, now pi (instead of py) dominates. 

driving frequencies (w < 2.5), E'max and i^min are given by the energies of the fourth 
instantaneous eigenstate, i.e. E'max = £'4(3/4) and £^min = E^^l/A), respectively, which 
means that the time-evolution is adiabatic. For higher driving frequencies, a series of 
resonances appear: i?max has sharp resonances at a; = 2.58, 3.45, 5.17, 6.83 and broad 
resonances with double peaks at a; = 10.19, 10.51 and u = 15.6, 16.4. On the other 
hand, -Emm shows considerably fewer resonances, sharp ones at a; = 3.97, 5.95 and a 
broad one at u = 11.89. The fact that i^min has considerably fewer resonances than 
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Figure 8. (Color online) Evolution of the energy E{t) for a driving frequency co = 100, 
when starting in the instantaneous ground state (first state) and in the fourth state. 
After roughly 60 driving periods, the energy E{t) saturates. 
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Figure 9. (Color online) Maximal and minimal energy reached when starting in the 
fourth state as a function of the driving frequency uj. i?max and -Emin show a series of 
resonances, which are related to a matching of the driving frequency uj to the mean 
energetic difference of instantaneous eigenstates, see tabled 

-E'max is immediately clear, since the only way that the energy, when starting in the 
fourth state (p4(0) = 1), can be lowered is to populate pi, whereas many different 
energetically higher lying instantaneous states can be populated, leading to an increase 

of -Emax- 

To gain an understanding of the above presented results, in particular of the 
resonances of -Emax and -Emin, we develop a Rabi-like model in the following section. 




max 
min 
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Table 2. (color online) Mean (averaged of the phase energy differences Ei^j 
together with Ei^^j/S and Ei^j/A for the first six states that share the same 

symmetry tt^ = iTy = +1. The energy differences which are related to the resonances of 
figure |9] between just two states are shaded in gray, the ones related to the large double 
peak at a; « 10.5 in yellow and the ones of the second double peak at w « 15 . . . w « 16 
in green. 
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5. Interpretation and Rabi-like model 

The aim is now to explain the resonances in figure [9] firstly by means of a population 
analysis of the instantaneous eigenstates and secondly via a Rabi-like model with time- 
periodic coupling. 

Since the instantaneous eigenstates and the corresponding energies in the elliptical 
billiard are quantized, only certain states can be populated when driving with a specific 
frequency u. In particular, we expect an efficient transfer between two states, when 
their energy difference exactly matches the driving frequency, i.e. when Ej — Ei = u. 
As we saw earlier, the energy Ei of an instantaneous eigenstate depends on the phase 
C, it is thus useful to define the mean energy difference -Ejoj between two instantaneous 
eigenstates as 

E,^, := (|E,(C) - i?,(C)|)c = \E^{Q - E,{Q\dC. (45) 

These mean differences together with Ei^j/2, Ei^j/?) and Ei^j/A are shown in table [2] 
for the first six states that share the same symmetry 11^^ = iVy = +1. Now if we start in 
the fourth state (p4(0) = 1), cf. figure [HI the only way to lower the energy is to populate 
the ground state (pi). Thus, when -E'io-4 ~ we expect a significant population of 



Resonant Population Transfer in the Time- Dependent Quantum Elliptical Billiard 19 

the ground state and consequently a decrease of -Emin- From table [2] i?i^4 = 11.92 and 
in figure [9], there is a broad resonance around u = 11.89. It is not surprising that the 
resonance is rather broad, since £'i.o-4 is just the mean energy difference, actually the 
difference between the first and the forth state varies between max{E4{() — Ei{()) = 14.6 
and min(£'4(C) — -Ei(C)) = 9.6, so transitions between the two states should be possible 
over a broad range around u = 11.89 and this is exactly what we observe. 

The transition between the first and the fourth state is a one photon transition in the 
sense that u = Ei^^, so one time the driving frequency matches the energy difference, 
yielding a broad resonance. On the other hand, there are also multi-photon transitions 
possible, such that n ■ u = -Ei^4, n = 2, 3, . . .. Since the efficiency of multi-photon 
processes drastically reduces with the number of photons involved, these higher-order 
resonances should be much sharper and are naturally harder to detect. The sharp 
resonances of -Emin at a; = 5.95 and u = 3.97 are the two and three photon resonances of 
-E'1^4, respectively, and the numbers are in excellent agreement with the corresponding 
values of table O 

In a similar fashion, the sharp resonances of -Emax can be explained: u = 
2.58, 3.45, 5.17 correspond to the four, three, and two photon transitions of E^^j, 
respectively and u = 6.83 to the three photon transition of E^^iq, again in excellent 
agreement with the results of table [2l Note that a population analysis of all the 
resonances confirms the so far presented considerations, e.g. at w = 6.83 besides p^, pio 
gets populated, whereas at a; = 5.17, besides ^4 now this is the case for pj. 

Applied to the broad double peaks of -Emax, a population analysis yields the 
following: at u = 10.19 a mixture of p^, pj and pio is present (also ^13, pi and P20 
are noticeably, however less than the first three ones populated), whereas at a; = 10.51 
mostly the states py, p4, Vn^ Pio and pi contribute. The yellow shielded fields 
of table |2] indicate that there is a multitude of one- to four-photon transitions in 
this frequency regime between exactly these states that according to the population 
analysis contribute to the resonance. Note that there are now also processes of the kind 
P4 — > — )■ pw possible, since the corresponding energy differences -^407 and -Ey^io 
are comparable. In contrast to the just described resonance, the broad double peak 
at a; = 15.6 and uo = 16.4 is mainly due to transitions between just two states. The 
corresponding population analysis yields p4 and pi3 for u = 15.6 and p^ and pis for 
u = 16.4, again in excellent agreement with the two-photon transition of -E'4^13 and the 
three-photon transition of £'4^18, see the greenish fields of table El 

For all the resonances, the population analysis additionally shows that the beating 
period Th increases dramatically when approaching a resonance, i.e. the population 
transfer from one instantaneous eigenstate to another becomes slower and slower the 
closer the driving frequency gets to the corresponding resonance. 

In the discussion of the resonances, we already used the projection onto the 
instantaneous eigenstates to analyze the results. With this picture in mind (see e.g. 
figure [31 where just two states are considerably populated), a very crude simplification 
is to interpret the driven elliptical billiard as a two-level system (two instantaneous 
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Figure 10. (Color online) Maximal and minimal energy reached when starting in 
fourth state as a function of the driving frequency uj now for the Rabi-like model. 
There is very good agreement with the full numerical simulations, see figure O 




eigenstates), where the two levels are coupled by a periodic driving with an unknown 
coupling strength a. Without the driving, the Hamiltonian of the unperturbed system 
in a suitable basis is then simply given by 

where Ei and E2 are the energies of the two eigenstates and \E'2; with time-evolution 
^O(t) = ^Og-*£„tM^ n = l,2. 

Next, we take into account the driving of the system, which leads to a periodic coupling 
of the two states. In the Hamiltonian, this leads to off-diagonal elements of the form 
asin{ujt). However, as can be seen in figure [H the energy itself of the instantaneous 
eigenstates depends on time. To this end, we introduce time-dependent energy shifts 
AEi{t) on the diagonal. In a good approximation, the AEi{t) are harmonic functions, 
which we obtain by fitting sine functions to the curves of figure [1] Overall, the model 
Hamiltonian is then 



n{t) = n° + 'H\t), 



where 



n'it) 



AEi{t) asm{ujt) 
asm{ujt) AE2(t) 



The standard ansatz for the time-dependent wave function "^{t) is a superposition with 
time-dependent coefficients of the unperturbed wave functions ^E'°(t), i.e. 

^{t) = c^{t)^Ut)+c2mi{t). 
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Plugging this ansatz into the time- dependent Schrodinger equation leads to a system of 
ordinary differential equations for the coefficients 

cJt) = \ \c2{t)asm{Lot)e'^^'-'^'^'/''' + Ci(t)AEi(t)l (46a) 
in 



C2{t) = [ci(t)asinMe'(^^-^^)*/'' + C2{t)AE2{t)] . (466) 



1 

ih 

Note that at this point, usually (e.g. in atomic systems) the rotating wave approximation 
is applied, however, in our case this would lead to incorrect results, so we have to solve 
the full system given by Eqs. (E]), which is only possible numerically. 

Not surprisingly, this simple 2-level model does not lead to correct results, e.g. 
the 2-level model cannot reproduce all the resonances of figure |9l The obvious reason 
is, as the population analysis shows, that already for driving frequencies in the range 
shown in figure |9] (a; < 17) more than just two states are involved. The straightforward 
generalization to n-levels yields: 

n{t) = n^ + H\t) 

where now V}{t) reads 

/ /\Ei{t) asinwt . . . \ 

a sin cut AE2(t) 
■H\t) = . ^ 

\ AE„(t) / 

or in other words 'Hjiit) = AEi{t) and for the off-diagonal elements we have 'H]j{t) = 
a sm{ujt). It is reasonable to assume that the coupling strength a is equal between all the 
levels, leaving us with just one free parameter in the model. Again, the time-evolution 
of the unperturbed states is given by 

and the superposition ansatz reads 

^{t) = j2^nm:{t), 

n 

yielding a set of coupled linear differential equations for the coefficients, which has to 
be solved numerically: 



-Ej)t/h 

ih 

j 



With this Rabi-like model, we calculate the time evolution of "^{t), with the same initial 
conditions as in the discussion of the resonances, i.e. we start in the fourth state, 
04(0) = 1. We take into account the first six states that can couple, i.e. the states with 
symmetries tTx = T^y = +1? namely 1, 4, 7, 10, 13 and 18, so we utilize a six-level Rabi- 
like model. Just like for the full numerical calculations, we compute -Emax and -E'min as 
a function of the applied driving frequency u. The result is shown in figure [TOl Overall, 
there is very good agreement with the results of the full simulations, cf. figure M In 
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particular, all resonances of figure IH] are reproduced. The additional resonances that 
appear in the Rabi-like model are four-, five- and even six-photon processes: u = 2.07 
corresponds to the five-photon transition of E^^t, u = 7.81 to the four-photon transition 
of -^4013 and u = 12.21 to the four-photon transition of -E'4<->i8- The resonance at 
u = 5.17 is superimposed with a very sharp resonance at a; = 5.20. This is the six- 
photon transition of £'4^13, additionally, we suppose that this resonance is enhanced 
by the combination of the two-photon transition E^^^^ with the four-photon transition 
E'r^is, again cf. table El 

The additional resonances are extremely sharp and this is one reason why they do 
not appear in figure M A much finer frequency scan would be necessary in the case of 
the full numerical simulations. Secondly, as already mentioned, the population transfer 
between two states becomes extremely slow when approaching a resonance (the beating 
period Tf, increases), i.e. the simulations have to be carried out over longer times in 
order to detect the very sharp resonances. 

The double peak at a; = 10.19, 10.51 of figure [9] appears in the Rabi-like model as 
a single, broad resonance, see figure [TOl The reason is that the time-dependent energy 
shifts AEi{t) of the model are obtained by fitting sine functions to the energy levels Ei{() 
of figure [U This fitting procedure does preserve the mean distances between the energy 
levels, in particular it reproduces the correct Ei^j, however, the minima of AEi{t) are 
lowered and the maxima of AEi[t) are raised compared to Ei{(). As a consequence, the 
resonances are broadened in the Rabi-like model, such that the mentioned double-peak 
merges, so they cannot be resolved within the model. Note that when fitting the Ei{() 
in a different fashion, such that the minima and maxima are preserved, the double peak 
does appear in the model, however this alternative fitting procedure changes the Ei^j, 
so all the resonances are shifted. 

The fact that the beating period Tf, increases drastically when approaching a 
resonance is also reproduced by the Rabi-like model, however, unlike for the appearance 
of the resonances itself, we do not have an intuitive explanation for these beatings. 
Nevertheless, the success of the Rabi-like model shows that we have identified the 
essential ingredients to explain the time-evolution of a wave function, when starting 
in an eigenstate, in the driven elliptical billiard. 

6. Conclusion and Outlook 

In this work, we have developed a numerical method to propagate an arbitrary 
initial state in the time-dependent elliptical billiard characterized by time-dependent 
boundary conditions. The method is based on a series of transformations of the 
original Hamiltonian, removing the time- dependent boundary conditions. Subsequently 
expanding the wave function in a suitable basis, we obtain a large system of coupled 
ordinary differential equations for the time-dependent expansion coefficients which can 
be solved by standard techniques. We then propagate eigenstates of the static system 
and investigate the time-evolution of their energy in the driven billiard. Unlike in 
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the corresponding classical system, the energy does not grow unboundedly but rather 
saturates, even for large driving frequencies, i.e. there is no Fermi acceleration in the 
quantum billiard. The maximal and minimal energy that is reached as a function of the 
applied driving frequency show a series of resonances. We analyze these resonances by 
projecting the wave function onto the instantaneous eigenstates, yielding the spectral 
composition of the states. At these resonances, population is periodically, with a 
period T^, transferred between the initial state and another one, either increasing or 
decreasing the energy of the wave function. This analysis shows that the resonances 
appear exactly when the driving frequency u matches the mean energetic difference 
between the involved instantaneous states (one photon transition) or multiples of the 
driving frequency (ra-photon transitions, which result in very sharp resonances). The 
beating period Tf, of the population transfer is much larger than the period T = 27t/uj 
of the driving and increases dramatically when approaching a resonance. To confirm 
our intuitive picture describing the appearance of the resonances, we develop a few- 
level Rabi-like model with time-periodic couplings, and time-dependent energy shifts 
on the diagonal, resembling the phase dependence of the energies of the instantaneous 
eigenstates. This model reproduces nicely the essential results of our full numerical 
simulations, in particular the mentioned resonances and the occurrence of the beating 
period Tj,. 

So far, we investigated the breathing mode of the elliptical billiard, i.e. the time- 
dependence of the semi-major and semi-minor axes is given by a{t) = + csin(a;t) 
and b(t) = Bq + csin(ci;t), respectively. A promising alternative would be to use a time- 
law obeying a{t)b{t) = const., since this preserves the area of the ellipse. As a result, 
the instantaneous energy levels exhibit a multitude of crossings, naturally posing the 
question whether a{t) could be adjusted in such a way that the momentum diffusion 
does not stop, yielding quantum Fermi acceleration in the driven elliptical billiard. 



Appendix A. Details of the ODE system 

The c^imn'l^) of ODE system, cf. Eq. (1371) . are given by 
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4 

where a{t) and b{t) are determined by the driving law of the eUiptical biUiard and the 
^nmn' defined as the following integrals over products of Bessel functions (for details 
of the derivation see EEI), where we use 



Inmn'{q,P,s) := i Jm{km,nr) Jm+p{km+q,n'r)r'' dv (A.4) 
^0 

and in the following we will omit the indeces nmn' for readability: 

Li = /(-2,0, 1) = 7(0,0,1) 

= 7(2,0,1) L4 = 7(0,2,1) 

L5 = 7(2,2,1) 7.6 = 7(4,2,1) 

L^ = 7(-4,-2,l) L« = 7(-2,-2,l) 

L*^ = 7(0, -2, 1) = 7(-l,0,0) 

= 7(1,0,0) L^^ = 7(-3,-2,0) 

7.13 = 7(_i, -2,0) 7.1^ = 7(1,2,0) 

7.15 = 7(3,3,0) Lis = 7(0, 0,3) 

7.1^ = 7(2,2,3) 7.1*^ = 7(-2,-2,3) 

The <iimn'(^) '^^^ t)e decomposed into time-independent parts f^mn' ^^"^ simple time- 
dependent functions gi{t) in the following way: 

= 9^'\t)f!rL + 9^'\t)fiZ.' (A.5) 
dlAt) = 9^'\t)f±, + g^'\t)f^i^, (A.6) 
= 9^'Ht)fi'l, + ^(^)(t)/il, (A.7) 



with 



9^'\t) = - (A.8) 
9 



^ ^ 

fia^t) fib^it) 

(2)(t) = (^a{t)d{t) + 6(t)6(t)) (A.9) 

= _ (A.IO) 

^?(^)(t)=/if(a(t)a(t)-Kt)Ki)) (A.ll) 



and 



J nmn' ' 



km,n'i-^nmn' '^^'nmn' ^ ^'nmn') 2/Cm^„/(m 1)7.1^^, 



8<-^m+l {kni,n) Jm+l {k-n'ni) 



8 Jm+l {km,n) Jm+l {kn' 
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r 16 

r(2) nmn' 

J nmn' '~ n J / 1, \ j J U 

^^m.+l \'^m,nj ^m+1 y^n 

f(-^) ._ 
J nmn' ' 



(3) km-2,n'{Lnmn' 2I/^^„, + ) + 6/Cm-2,n' (^^ l)-^nmn' 



1.Q Jm+li^kjjin) .J^—l {km—2,n' 

2km-2,n'{m-3)Lll^^, 



lQJm-\-l(^k^n) J m—1 {km—2,r 
r 18 

M) ^nmn' 

J nmn' ' 



4t/m+l {km,n) Jm—1 {km—2,n 

J nmn' ' 



(5) ^m+2,n'i'^nm.n' '^'^nmn' + ^^nmn') + 6/Cm+2,n'('^ + ^)'^nm.n' 

16 Jm+1 {km,n) Jm+d,{km+2,n' ) 
'2km+2,n'ijn + 3)L]^rnn' 



1 6 J m+1 ( km,n) J m+3 ( ^m+2,n' ) 



r 17 

._ nmn' j^2^ 

nmn ^Jm.+l{km^n) Jm+z{km+2,n') 
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